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The  finned  tube  contactor  of  a  liquid  desiccant  (triethylene 
glycol)  dehumidifier  system  is  modeled  with  a  fifth  order  differential 
equation  to  simulate  the  simultaneous  heat  and  mass  transfer  associated 
with  conditioning  the  air.   Potential  energy  saving  aspects  of  the  sys- 
tem are  discussed.   Computer  simulation  results  presented  on  psychro- 
metric  charts  show  that  significant  performance  is  possible  with  contac- 
tors of  one  half  the  recommended  size.   The  cycle  may  be  regenerated  at 
120  to  140°F  depending  on  the  system  design  requirements. 

Figure  4  may  be  used  to  estimate  system  performance. 
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NOMENCLATURE  AND  UNITS 


Abbreviations 


LDD  Liquid  Desiccant  Dehumidif ication 

TEG  Triethylene  Glycol 

English  Letter  Symbols 

a  Cooling  area  per  unit  volume  for  heat  transfer, 

ft2/ft3 

a^  Interfacial  area  per  volume  for  heat  transfer, 

ft2/ft3 

a^  Interfacial  area  per  volume  for  mass  transfer, 

ft2/ft3 

c,C  Heat  capacity  at  constant  pressure,  BTU/lb 

Cq  Humid  heat,  BTU  (of  water  vapor-gas  mixture)/ 

(lb  dry  gas) (°F)  defined  by  Eq.  (A. 22) 
C  Humid  heat,  BTU  (of  TEG-water  solution) /(lb  TEG) 

(°F)  defined  by  Eq.  (A. 33) 

D  Root  diameter  of  tube,  ft 

r 

D  Diffusion  coefficient  for  water  vapor  in  air, 


defined  by  Eq.  (C.l) 

2 
Acceleration  due  to  gravity,  ft/hr 

Superficial  gas  (dry  air  plus  water  vapor)  mass 

velocity,  lb/(hr)(ft2) 


Gas  mass  velocity  at  minimum  cross  section, 


Superficial  solvent  (dry  air)  mass  velocity, 


h  Contact  coefficient  of  heat  transfer  between  gas 

and  liquid,  BTU/(hr) (ft2) (°F)  defined  by  Eq.  (C.2) 

h'  Contact  coefficient  of  heat  transfer  between  gas 

and  liquid  corrected  for  simultaneous  mass  trans- 
fer BTU/(hr)(ft2)(°F)  defined  by  Eq.  (A. 15)  and 
(A. 16) 

h  Contact  coefficient  of  heat  transfer  between 

liquid  and  cooling  fin,  BTU/(hr)(ft  )(°F)  defined 
by  Eq.  (C.13) 

h  Conduction  coefficient  of  heat  transfer  in  liquid, 

BTu7(hr)(ft2)(°F)  defined  by  Eq.  (C.ll) 

H  Enthalpy,  BTU/lb 

Hp  Heat  of  dilution  (mixing),  BTU/lb 

k  Thermal  conductivity,  BTU/(hr) (f t) (°F) 

1  Effective  fin  length,  in. 

L  Superficial  liquid  (TEG  plus  water)  mass  velocity, 


Superficial  solvent  (TEG)  mass  velocity, 


defined  by  Eq.  (A. 14) 
p  *  Partial  pressure  in  equilibrium  with  a  liquid 

solution,  mm  Hg 
p  Total  pressure,  mm  Hg 

P  Fin  pitch,  fins/in;  Chapter  II:   pressure 

Q  Heat  removed  by  cooling  coil,  BTU/hr  defined  by 

Eq.  (A. 28) 
Q  Exothermic  heat  of  dilution,  BTU/hr  defined  by 

Eq.  (A. 30) 

q  Heat  transfer  flux  for  sensible  heat,  BTU/(hr) 

s 

(ft2) 

R  Chapter  II:   universal  gas  constant; 

Appendix  B:   bond  graph  resistance 

s  Spacing  between  fins,  in. 

2 
S  Interfacial  surface,  ft  /f 

t  Thickness  of  aluminum  cooling  fin,  in;  tempera- 

ture, °F 
T  Chapter  II ,  temperature ,  °F 

T  Horizontal  tube  spacing,  in. 

T  Vertical  tube  spacing,  in. 

sv  r    ©» 

U  Over-all  coefficient  of  heat  transfer  between  the 

2 
liquid  and  cooling  water,  BTU/(hr)(ft  ) (°F) 

V  Water  velocity  through  tube,  ft/sec 
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x  Mass  fraction  concentration  in  the  bulk  liquid, 

lb  water/ (lb  water  +  lb  TEG) 

X  Mass  ratio  concentration  in  the  bulk  liquid, 

lb  water/lb  TEG 

y  Mass  fraction  concentration  in  the  bulk  gas, 

lb  water  vapor/ (lb  water  vapor  +  lb  dry  air) 

y*  Mole  fraction  concentration  in  the  gas  in  equilib- 

rium with  the  bulk  liquid,  moles  water  vapor/ 
(moles  water  vapor  +  moles  dry  air) 

y  Mole  fraction  concentration  in  the  bulk  gas, 

m 

moles  water  vapor/ (moles  water  vapor  +  moles  dry 

air) 
Y  Mass  ratio  concentration  in  the  bulk  gas,  lb 

water  vapor/ (lb  dry  air) 
Z  Height  of  finned  tube  contactor  (parallel  to 

flow  direction),  ft. 

Greek  Letter  Symbols 

T  Liquid  loading  per  unit  width  of  fin,  lb/(hr)(ft) 

6  TEG  solution  thickness,  ft. 

n.  Fin  efficiency,  dimensionless 

X  Latent  heat  of  vaporization,  BTU/lb 

u  Viscosity  at  bulk  temperature,  lb/(hr)(ft) 

3 
p  Density,  lb/ft 

o  Ratio  of  free  flow  area  to  face  area,  dimension- 

less 


V 


Chapter  II:   specific  volume,  ft  /lb 


Dlmensionless  Groups 

Pr 

Sc 

Subscripts 

A 

AL 

B 

C 

D 

f 

i 

1,L 

o 

Bond  Graph  Symbols 
C 


Q 
R 

Ra 
Ri 


Prandtl  number  (uc  /k) 
P 

Schmidt  number  (u/p  D  ) 

Substance  A,  water  vapor 

Substance  A,  liquid  water 

Substance  B,  dry  air 

Substance  C,  TEG 

Substance  D,  water  plus  TEG  solution;  dilution 

Cooling  water 

Gas 

Interface 

Liquid 

At  the  reference  temperature 

Capacitance 

Mass  capacitance 

Thermal  capacitance 

Mass  transfer  flux 

Heat  flux 

Resistance 

Diffusion  resistance 

Interfacial  resistance  to  mass  transfer 


R  Reaction  resistance 

r 

S^,  Heat  of  dilution  flow  source 
fd 

Sf  Latent  heat  flow  source 

t  Temperature 

T  Transformer 

0  Parallel  junction,  flows  sum  to  zero 

1  Series  junction,  efforts  sum  to  zero 
Power  bond 

> Signal  bond 


I.   INTRODUCTION 

Definitions 

Liquid  desiccant  dehumidif ication  is  the  reduction  of  the  water 
vapor  content  of  a  given  volume  of  air  by  a  liquid  sorbent  which  has  the 
property  of  extracting  and  holding  water  vapor  brought  into  contact  with 
it.   The  terms  sorbent  and  desiccant  refer  to  materials  with  a  large  ca- 
pacity for  moisture  compared  to  their  volume  and  weight.   The  amount  of 
water  held  by  a  desiccant  will  increase  or  decrease  depending  on  whether 
the  vapor  pressure  of  water  above  the  desiccant  is  less  than  or  greater 
than  the  partial  pressure  of  water  vapor  in  the  air  in  contact  with  it. 

Applications 

Dehumidif ication,  whether  by  liquid  desiccants  or  by  other  means, 
has  important  commercial  applications.   These  include  preserving  ships 
and  other  surplus  equipment  which  would  otherwise  deteriorate,  and  pre- 
venting condensation  damage  to  cargo  and  internal  structures  of  ships' 
holds  which  transport  dry  cargo. 

Liquid  desiccant  dehumidification  (LDD)  was  investigated  in  the 
1950' s  and  found  to  be  operational  [7], [22].   Currently  LDD  equipment 
designed  for  low  pressure  steam  regeneration  is  commercially  available 
and  is  used  in  the  drying  of  natural  gas  prior  to  pipeline  transporta- 
tion.  Both  hospitals  and  the  food  processing  industry  use  LDD  not  only 
for  temperature  and  humidity  control  but  also  for  the  bactericidal  prop- 
erties of  the  family  of  glycol  desiccants.   Accepted  high  equipment  cost 
prevents  wider  application  of  LDD  for  the  general  purpose  of  space  de- 
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humidification  or  space  cooling.   These  costs  are  on  the  order  of  $10 
per  square  foot  of  contactor  in  typical  LDD  systems  which  remove  about 
500  pounds  of  water  per  hour  with  dehumidification  and  regeneration  con- 
tactor areas  of  5000  sq.  ft.  each  [27]. 

Potential  Savings 

If  the  potential  to  regenerate  an  LDD  cycle  at  low  temperatures  is 
exploited  and  waste  heat  or  solar  energy  is  utilized,  then  potential 
energy  savings  can  be  applied  against  acquisition  cost.   Such  cost  may 
be  minimized  by  rational  rather  than  empirical  equipment  selection. 

A  review  of  current  space  dehumidification  and  cooling  practice 
will  indicate  the  potential  energy  savings  of  LDD. 

In  conventional  air  conditioning  systems  the  air  is  cooled  signifi- 
cantly below  the  desired  temperature  to  effect  moisture  removal.  (Refer 
to  Figure  1). 


CHILLED  WATER 

ENERGY 

1 

WARM.  MOIST  AIR 

COOLING 
COIL 

VERY  COOL,  DRY  AIR 

REHEAT 
COIL 

COOL  DRY  AIR 

T 

Figure  1.   Conventional  Air  Conditioner  with  Reheat 
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Additional  energy  is  used  as  this  very  cool,  dry  air  is  then  reheated 
to  the  desired  space  temperature.   The  popularity  of  this  system  stems 
from  the  previously  low  cost  of  energy  and  the  simple  control  system  re- 
quired.  One  dependent  variable,  outlet  humidity  ratio,  is  uniquely  con- 
trolled by  the  chilled  water  temperature  through  the  condensation  of  ex- 
cess moisture  above  the  dew  point  temperature.   The  other  dependent  vari- 
able, desired  space  temperature,  is  now  uniquely  controlled  at  a  constant 
humidity  by  the  quantity  of  reheat  energy  added. 

The  control  system  complexity  increases  when  an  attempt  is  made  to 
save  reheat  energy  by  mixing  the  very  cool  dry  air  with  warm  moist  air 
bypassed  around  the  cooling  coil.   (Refer  to  Figure  2.) 


CHILLED  WATER 


WARM,  MOIST  AIR 


COOLING 
COIL 


BYPASS  AIR 


Figure  2.   Conventional  Air  Conditioner  with  Bypass  Mixing 

Both  of  the  two  independent  variables,  chilled  water  temperature  and  bypass 
air  flow  rate,  simultaneously  modify  both  of  the  dependent  variables, 
outlet  air  temperature  and  humidity  ratio. 

The  LDD  system  can  save  the  energy  of  excess  chilling  for  moisture 
removal  through  its  ability  to  dehumidify  at  ambient  temperatures.  Re- 
moval of  the  latent  heat  of  condensation  is  the  primary  function  of  the 
cooling  water.   (Refer  to  Figure  3.) 


CONCENTRATED 
LIQUID  DESICCANT 


COOLING  WATER 


WARM  MOIST  AIR        AIR-LIQUID      COOL  DRY  AIR 
*    CONTACTOR  ' 


DILUTE 

LIQUID  DESICCANT 


Figure  3.   Liquid  Desiccant  Dehumidif ication  System 


The  dehumidif ication  contactor  of  the  LDD  system,  a  three  fluid  heat  ex- 
changer in  which  mass  transfer  between  two  of  the  fluids  occurs,  can 
supply  cool,  dry  air  without  further  heating  or  cooling.   Thus,  the  re- 
heat energy  can  also  be  saved. 

In  the  LDD  contactor  the  two  dependent  variables,  air  temperature 
and  humidity  ratio,  are  controlled  by  the  independent  variables,  liquid 
desiccant  concentration  and  cooling  water  temperature.   A  tradeoff  for 
the  potential  energy  savings  of  excess  chilling  and  unnecessary  reheat 
is  a  control  system  of  sophistication  similar  to  that  of  the  bypass  sys- 
tem shown  in  Figure  2. 

The  LDD  cycle  compares  favorably  to  an  alternative  technique  for 
space  cooling  such  as  the  lithium  bromide  absorption  cycle.   While  the 
absorption  cycle  is  closed  and  the  desiccant  cycle  is  open,  both  systems 
have  coefficients  of  performance  on  the  order  of  0.6  [11].   The  absorp- 
tion cycle  requires  regeneration  temperatures  around  200°F,  while  the 
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LDD  cycle  can  be  regenerated  at  temperatures  around  140 °F.   Of  these 
two  systems  the  absorption  system  has  been  the  subject  of  considerable 
analysis  and  simulations  [9];  however,  there  is  no  well-established  tech- 
nique for  space  cooling  or  dehumidification  with  low  temperature  regener- 
ation. 

Purpose 

In  view  of  the  above,  a  model  for  analysis  of  the  liquid  desiccant 
dehumidifier  contactor  of  Figure  3  is  presented  as  the  basis  for  a  ra- 
tional design.  The  primary  function  of  the  model  is  to  assure  that  the 
contactors  are  not  over-specified  so  that  acquisition  costs  may  be  kept 
as  low  as  possible.  The  contactors  (dehumidification  and  regeneration) 
represent  a  significant  portion  of  the  total  cost  of  an  LDD  system. 

Empirical  equipment  selection  of  the  contactors  is  standard  prac- 
tice and  can  result  in  overspecification  to  ensure  proper  operation. 
Manufacturers'  application  recommendations  specify  equipment  which  oper- 
ates with  the  exit  conditions  of  air  and  liquid  temperatures,  humidity 
ratio,  and  desiccant  concentration  at  equilibrium. 

By  contrast,  rational  design  implies  more  than  black  box  modeling 
of  the  contractor.   The  computer  model  of  the  contactor  which  was  devel- 
oped is  based  on  a  complete  set  of  energy  balances,  material  balances, 
rate  equations,  and  equilibrium  relationships.   Through  computer  simu- 
lation of  the  heat  and  mass  transfer  in  the  contactor,  the  rate  at  which 
the  temperatures,  humidity  ratio,  and  desiccant  concentration  approach 
equilibrium  as  a  function  of  flow  distance  can  be  studied.   With  a  speci- 
fic design  the  contactor  flow  length  and  associated  cost  can  be  compared 


18 

to  the  incremental  heat  and  mass  transfer  of  additional  contactor  flow 
length  increments. 
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II .   PSYCHROMETRY 

The  Psychrometric  Chart 

Here  psychrometry  is  concerned  with  the  determination  of  the  prop- 
erties of  the  mixture  of  air  and  water  vapor  as  the  mixture  flows  through 
the  contactor.   These  properties  are  conveniently  diagrammed  in  an  abbre- 
viated psychrometric  chart,  Figure  4,  which  plots  dry  bulb  temperature, 
humidity  ratio,  wet  bulb  temperature,  saturation  temperature,  and  rela- 
tive humidity  for  the  air-water  vapor  mixture. 

Dry  bulb  temperature  represents  the  air  and  liquid  contact  tempera- 
ture and  is  the  abscissa. 

Humidity  ratio  (H.R.)  is  the  pounds  of  water  carried  by  one  pound 
of  dry  air.   Alternatively,  grains  per  pound  of  dry  air  (grains)  is  used 
where  7000  grains  equal  one  pound  of  water.   The  ordinate  shows  both 
humidity  ratio  and  grains. 

Relative  humidity  is  the  ratio  of  partial  pressure  of  the  water 
vapor  to  the  vapor  pressure  of  water  at  the  air  temperature  and  is  indi- 
cated by  the  dotted  lines.   Relative  humidity  is  expressed  in  percent. 

Wet  bulb  temperature  lines  run  upward  and  to  the  left  to  the  satu- 
ration curve  and  very  nearly  represent  the  path  followed  by  adiabatic 
(spray)  cooling. 

Dew  point  temperature  is  the  temperature  to  which  the  mixture  must 
be  cooled  to  become  saturated  and  is  found  by  following  a  horizontal  line 
to  the  saturation  temperature  curve. 

Saturation  temperature  represents  the  maximum  moisture  the  air  can 
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hold  before  fogging  ensues. 

Enthalpy  and  specific  volumes  are  plotted  on  more  detailed  charts 
[26].   Any  two  of  these  parameters  specify  the  air-water  vapor  mixture 
and  allow  all  the  others  to  be  determined.   The  remaining  parameter, 
humectant  value,  is  discussed  later  under  its  own  section. 

Dry  bulb  temperature  and  humidity  ratio  were  chosen  as  the  two  de- 
pendent variables  of  interest  because:   (1)  dry  bulb  temperature  is 
easily  measured  and  has  significance  to  people;  and  (2)  the  humidity 
ratio  directly  represents  the  water  content  of  the  air  and  is  based  on 
one  pound  of  dry  air  which  is  a  constant  in  any  psychrometric  process. 

Psychrometric  Chart  Processes 

The  process  of  heating  or  cooling  moist  air  above  the  dew  point  is 
represented  by  a  horizontal  line  since  the  humidity  ratio  is  unchanged. 
This  process  is  known  as  sensible  heating  or  cooling. 

The  process  of  cooling  below  the  dew  point  entails  condensation  of 
vapor  and  is  represented  by  a  horizontal  line  extending  left  to  the  satu- 
ration curve  where  the  process  follows  the  curve  to  a  lower  humidity 
ratio.   This  is  roughly  the  dehumidif ication  process  of  conventional  air 
conditioning. 

A  downward  vertical  line  on  the  chart  represents  an  isothermal  dehu- 
midif ication  process  or  latent  cooling.   This  may  be  accomplished  directly 
by  the  LDD  system  with  cooling  to  remove  only  the  latent  heat  of  conden- 
sation. 

Liquid  Desiccants 

A  liquid  desiccant  removes  water  vapor  from  the  air  at  a  rate  direc- 
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tly  proportional  to  the  difference  between  the  partial  pressure  of  the 
water  vapor  immediately  above  the  liquid  desiccant.   The  partial  pres- 
sure of  water  vapor  in  the  air  is  reduced  and  the  liquid  desiccant  is 
diluted  in  the  process  until  equilibrium  occurs.   The  process  may  be  re- 
versed by  heating  the  desiccant.   That  is,  the  liquid  desiccant  is  con- 
centrated as  the  partial  pressure  of  water  vapor  in  the  air  increases. 
The  vapor  pressure  of  the  liquid  desiccant  is  a  function  of  temperature 
and  concentration. 

Humectant  Values 

Equilibrium  humectant  values  are  the  percent  by  weight  of  a  liquid 
desiccant  in  water  solution  required  to  maintain  vapor  pressure  equili- 
brium in  contact  with  air  of  various  temperatures  and  humidity  ratios. 
The  humectant  values  for  triethylene  glycol  (TEG)  were  added  to  the 
standard  psychrometric  chart  to  produce  a  convenient  diagram  revealing 
the  relationship  between  TEG  concentration  and  the  air-water  vapor  mix- 
ture at  equilibrium.   The  humectant  values  on  Figure  A  were  developed 
from  graphs  of  the  vapor  pressure  of  water  over  triethylene  glycol  of 
various  weight  concentrations  and  temperatures  [28]  [37],   These  equilib- 
rium humectant  values  give  an  indication  of  the  TEG  concentration  necess- 
ary to  dehumidify  to  a  specified  outlet  condition  with  a  constant  TEG 
concentration. 

The  solid  lines  which  are  the  TEG  equilibrium  humectant  values  rep- 
resent slightly  different  TEG  weight  concentrations  in  the  90%  to  99% 
range  depending  on  the  data  base:   Union  Carbide  [37]  or  Dow  Chemical  [281 
A  third  source  [12],  presented  data  for  only  three  concentrations,  40%, 
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60%,  and  80%  and  these  were  consistent  with  the  first  two  sources  within 
5%.   The  vapor  pressures  presented  by  Union  Carbide  and  Dow  Chemical 
differ  by  approximately  20%  at  90%  concentration  and  100%  at  98%  concen- 
tration. Union  Carbide  data  presented  the  higher  vapor  pressures.   The 
more  conservative  humectant  values  during  dehumidification  are  derived 
from  the  Union  Carbide  data  while  the  more  conservative  humectant  values 
during  regeneration  are  derived  from  the  Dow  Chemical  data. 

In  the  area  of  Figure  4  in  which  dehumidification  takes  place,  70 
to  80  °F  and  0.010  H.R.,  the  change  in  TEG  humectant  value  with  respect 
to  temperature  is  about  1%/°F.  At  a  0.025  H.R.  and  120  to  130  °F  where 
regeneration  can  take  place,  the  ratio  is  about  0.2% /°F  which  indicates 
that  TEG  regeneration  should  not  be  particularly  sensitive  to  regener- 
ation temperatures. 

Vapor  Pressure  Prediction 

An  equation  for  the  equilibrium  partial  pressure  above  the  liquid 
desiccant  is  necessary  for  the  computer  simulation.   In  addition  to 
accuracy,  a  suitable  equation  will  have  continuous  derivatives  over  its 
range  and  will  be  relatively  efficient  with  respect  to  calculation  time. 
Two  relationships  with  these  characteristics  which  predict  vapor  pres- 
sure as  a  function  of  TEG  contact  temperature  and  TEG  weight  concen- 
tration were  developed  with  the  aid  of  a  curve  fitting  program  and  a 
programmable  calculator: 

From  Dow  Chemical  data, 
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From  Union  Carbide  data 
p*  [mm  Hg]  =  exp 
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(2.2) 


0.26921 
where  T  =  contact  temperature, 


Both  relationships  are  accurate  within  about  3%  within  the  ranges 

40  <  T  <  300,  75  <  C  <  99  and  are  slightly  less  accurate  for  0  <  C  <   75, 

Eq.  (2.2)  was  used  exclusively  in  the  simulation. 

Vapor  Pressure  Correlation 

The  method  used  in  synthesizing  equations  (2.1)  and  (2.2)  is  based 
on  the  Clausius-Clapeyron  equation  which  relates  the  slope  of  the  vapor- 
pressure  curve  to  the  latent  heat  of  vaporization 

dP  _     X  (2.3) 

dT   T(VG  -  VL) 

where  v  =  specific  volume  of  the  saturated  vapor 


L 


specific  volume  of  the  saturated  liquid 


X  =  latent  heat  of  vaporization 
Using  the  ideal-gas  law  and  neglecting  V  in  comparison  with  V 

dtaP-ffi-i«  (2.4) 

RT 

If  X  is  fairly  constant  over  the  temperature  range 

In  P  =  -  —  +  constant  (2.5) 
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Equation  (2.5)  suggests  that  a  plot  of  log  P  against  —will  be  linear 
over  small  temperature  ranges.   Antoine's  equation  published  in  1888 
allows  for  slight  variation  of  A  with  temperature  [26] 

(2.6) 

and  is  probably  the  best  vapor-pressure  correlation  equation  for  a  simple, 
three-constant  (a,  b,  and  c)  equation.   Equation  (2.6)  is  the  basis  for 
equations  (2.1)  and  (2.2). 


III.   LIQUID  DESICCANT  DEHUMIDIFICATION  SYSTEM 

Operating  Cycle 

One  system  for  liquid  desiccant  dehumidif ication   (LDD)   is  shown 
in  Figure  5.   The  dehumidif ier  section  on  the  left  accepts  warm  moist 
air  from  the  atmosphere  and  contacts  the  air  with  the  liquid  desiccant 
in  the  dehumidif ier  contactor.   Intimate  contact  between  the  air  and 
desiccant  over  a  large  area  occurs  through  the  formation  of  a  thin  film 
of  desiccant  as  it  is  sprayed  over  the  extended  surface  contactor.   The 
heat  of  condensation  and  the  heat  of  dilution  flow  through  the  liquid  to 
the  cooling  water  which  can  be  sea  water,  well  water,  or  water  from  a 
cooling  tower.   The  quantity  of  cooling  water  required  is  strongly  de- 
pendent on  the  sensible  and  latent  heat  removed  from  the  air  by  the 
liquid  desiccant  and  weakly  dependent  on  the  desiccant  circulation  rate 
between  the  dehumidif ier  and  regenerator. 

The  air  leaves  the  contactor  flowing  downward  with  entrained  desic- 
cant and  loses  most  of  the  liquid  to  the  sump  while  crossing  it.   The  re- 
mainder of  the  entrained  liquid  is  eliminated  as  the  cool,  dry  air  flows 
upward  to  the  controlled  space. 

The  dry  bulb  temperature  of  this  cool,  dry  air  is  weakly  dependent 
on  the  water  removed  from  the  air  but  is  strongly  dependent  on  the  temp- 
erature of  the  liquid  desiccant.   Accordingly,  the  air  flowing  through 
the  dehumidifier  can  be  warmed  or  cooled  as  it  is  dried  depending  on  the 
temperature  of  the  liquid  desiccant  and  the  temperature  of  the  cooling 
water. 

Concentration  of  the  liquid  desiccant  is  more  easily  achieved  with 
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dry  air  which  would  normally  be  exhausted  from  the  dehumidified  space, 
but  can  also  be  achieved  with  outside  air.   The  heat  of  vaporization 
and  dilution  is  supplied  through  the  heating  coil  in  the  regenerator 
contactor  to  the  liquid  desiccant  as  it  contacts  the  relatively  dry  air. 
Heating  evaporates  the  liquid  desiccant  though  at  a  much  reduced  rate 
than  that  of  the  water.  A  reflux  coil  maintained  at  about  90  °F  is  used 
as  a  condenser  for  the  desiccant  to  prevent  excessive  loss  in  the  hot 
exhaust. 

Figure  5  also  shows  the  four  pumps  necessary  for  the  liquid  desic- 
cant movement — two  spray  pumps  and  two  feed  pumps.   The  liquid-to-liquid 
heat  exchanger  between  the  conditioner  sump  which  is  slightly  warmer 
than  the  cooling  water,  and  the  concentrator  sump  which  is  slightly 
cooler  than  the  heating  water,  improves  the  system  thermal  efficiency. 
An  air-to-air  heat  exchanger  between  the  hot  wet  exhaust  and  the  cooler 
inlet  air  of  the  regenerator  will  also  increase  the  system  thermal  effi- 
ciency. 

Figure  6  shows  a  similar  system  with  a  combined  sump  which  requires 
only  one  pump.  While  this  system  is  more  compact  and  uses  fewer  pumps, 
thermal  efficiency  can  be  degraded  since  the  sump  operates  at  a  mean 
temperature.  With  equal  dehumldifier  and  regenerator  desiccant  recircu- 
lation rates  the  sump  operates  as  a  counter  flow  heat  exchanger  with  an 
effectiveness  of  0.5,  or  less  than  is  normally  expected  from  the  counter 
flow  heat  exchanger  of  Figure  5. 

In  either  system  the  solution  in  the  sump  will  require  treatment  as 
it  will  concentrate  solids  absorbed  from  the  air. 
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Air-Desiccant  Contactor 

The  contactor  is  the  part  of  the  system  in  which  the  air  is  brought 
in  contact  with  the  liquid  desiccant  for  dehumidif ication  or  regener- 
ation.  The  contactor  should  have  high  heat  and  mass  transfer  coeffici- 
ents but  yet  have  relatively  low  air  flow  resistance.   The  contactor 
should  also  be  light,  compact,  inexpensive  and  retain  its  heat  and  mass 
transfer  characteristics  after  extended  use. 

Although  the  contactor  may  consist  of  a  simple  spray  with  the  re- 
sulting air-desiccant  interface  being  the  surface  of  droplets,  entrap- 
ment would  be  excessive  since  the  finely  divided  liquid  is  difficult  to 
separate  from  the  air.   In  addition,  it  would  be  difficult  to  remove  the 
heat  of  condensation  and  dilution.   A  drop  spray  onto  extended  surface 
finned  tubes,  Figures  7  and  8,  prevents  atomization  of  the  liquid  desic- 
cant and  has  a  better  mass  transfer  capacity  per  volume  than  a  simple 
spray  [33]. 
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Figure  7.   Four  Row  Finned  Tube  Contactor  Dimensions 


Dimensions  for  Figure  7 
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Figure  8.   Dimensions  for  Determining  Free  Flow  Area  of  Contactor 
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The  cooling  water  flows  through  the  tubes  in  a  parallel-counter- 
flow  configuration  with  an  even  number  of  tube  passes,  while  the  air  and 
liquid  desiccant  flow  over  the  extended  surface.   Thus,  the  dehumidifi- 
cation  contact  surface  and  cooling  coil  are  conveniently  packaged  in  one 
unit.   The  cooling  water  serves  to:   (1)  keep  the  desiccant  cool  and  thus 
lower  the  vapor  pressure;  and  (2)  cool  the  air. 

The  regenerator  contactor  has  the  same  physical  characteristics 
and  performance  criteria  as  the  dehumidifier  contactor,  although  the  size 
may  differ  because  of  different  air  flow  and  liquid  desiccant  recircu- 
lation rates.   Therefore,  the  computer  model  presented  later  simulates 
the  heat  and  mass  transfer  of  both  the  dehumidifier  and  regenerator  con- 
tactors. 

The  air  flow  can  be  either  parallel  with  or  counter  to  the  liquid 
desiccant  flow  direction,  depending  on  the  particular  application.  Due 
to  the  dilution  of  the  liquid  desiccant  as  it  flows  through  the  contac- 
tor, the  exit  humidity  ratio  is  not  as  low  as  can  be  achieved  with  the 
counter  flow  arrangement.  However,  the  downward  parallel  flow  configu- 
ration does  not  inhibit  spraying  the  liquid  on  the  contactor  and  allows 
some  separation  of  the  liquid  and  air  by  a  reversal  of  flow  direction. 

The  parallel  flow  configuration  is  best  suited  for  the  processing  of 
large  volumes  of  air,  while  the  counter  flow  configuration  can  deliver 
the  very  dry  air  needed  in  manufacturing  operations  such  as  the  handling 
of  hygroscopic  materials.   The  parallel  flow  configuration  represents 
our  physical  system. 
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Other  System  Equipment 

Other  equipment  items,  such  as  the  main  circulating  blowers  or 
liquid  sorbent  pumps,  should  be  selected  on  the  basis  of  good  engineer- 
ing design  with  cognizance  of  the  difficulty  of  containing  TEG  because 
of  its  propensity  to  leak  past  defective  fittings  or  packings. 

LDP  System  Advantages 

1.  Low  temperature  regeneration  as  opposed  to  other  thermal  cycle  air 
conditioners  is  possible  and  allows  more  efficient  operation  of  solar 
collectors. 

2.  Regeneration  is  not  very  sensitive  to  heating  water  temperature,  un- 
like other  thermal  cycle  air  conditioners. 

3.  Solar  or  recovered  waste  power  can  be  used  to  regenerate  the  liquid 
desiccant . 

4.  The  system  is  simple  in  construction  and  operation  since  the  cycle 
is  without  high  pressure  liquids  or  gases. 

5.  Precise  control  of  air  temperature  and  humidity  results  from  a  speci- 
fic TEG  concentration  and  cooling  water  temperature. 

6.  Co-sorption  benefits  include  the  removal  of  bacteria,  odors  and 
various  airborne  irritants  [1]. 

7.  The  expense  of  a  separate  cooling  coil  can  be  saved  or  a  smaller 
after  cooler  can  be  specified  if  chilled  water  is  used  in  the  dehumidi- 
fication  contactor. 

8.  Excess  dehumidification  capacity  when  obtained  from  low  cost  waste 
heat  or  solar  energy  may  be  advantageously  converted  to  sensible  cooling 
by  spray  cooling.   The  spray  process  cools  and  humidifies  essentially 
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along  the  lines  of  constant  wet  bulb  temperature  in  Figure  4. 
9.   In  ocean  engineering  applications  dew  point  temperatures  on  the 
order  of  10  °F  below  the  structure  temperature  are  required  to  prevent 
condensation  [10].   The  combination  of  LDD  and  cold  sea  water  can  produce 
sub-freezing  dew  point  temperatures  as  the  antifreeze  qualities  of  TEG 
prevent  contactor  frost  accumulation. 


IV.   METHOD  OF  ANALYSIS 

Mathematical  Model 

The  analysis  of  the  contactor  is  based  on  the  development  of  a 
mathematical  model  suitable  for  a  computer  simulation.   The  substance 
of  the  model  developed  is  presented  here;  detailed  information  may  be 
found  in  the  appendices. 

The  model  of  the  contactor  is  a  fifth  order  differential  equation 
composed  of  three  first  order  differential  equations  for  heat  transfer 
and  two  first  order  differential  equations  for  mass  transfer.   The 
three  heat  transfer  state  variables  are  the  gas,  liquid,  and  cooling 
water  temperatures  as  shown  in  Figure  9. 
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Figure  9.   Locations  of  Representative  Temperatures  and  Sources 
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TEG  Solution  thickness  0.007  to  0.010  in 
Fin  thickness  0.010  in 

Tube  wall  thickness  0.035  in 

Tube  diameter  0.750  in  O.D. 

Dimensions  for  Figure  9. 

The  two  mass  transfer  state  variables  are  the  masses  of  the  gas,  (water 
vapor  plus  dry  air),  and  liquid,  (water  plus  TEG). 

All  differentials  are  taken  with  respect  to  flow  distance  in  a 
steady  state  model.   This  is  a  reasonable  assumption,  as  start  up  trans- 
ients will  be  negligible  within  a  minute  or  two  of  operation  from  a  de- 
termination of  the  system  RC  time  constants. 

Actual  variables  used  in  the  program  for  the  mass  state  variables 
are  the  superficial  mass  velocities  of  the  gas  and  liquid.   Thus,  the 
five  independent  variables  in  this  model  are  t  ,  t  ,  t  ,  G,  and  L.   All 
the  other  variables  are  dependent  variables  and  can  be  calculated  from 
known  relationships  with  knowledge  of  the  five  independent  variables. 
Two  distinct  methods  were  used  to  develop  the  mathematical  model. 
Appendix  A  presents  a  derivation  using  the  traditional  method  for  analyz- 
ing mass  transfer  operations  based  on  material  and  energy  balances  in 
control  volumes  of  a  differential  section  of  the  contactor.   Appendix  B 
uses  a  general  method.   It  develops  a  bond  graph  model  of  the  contactor 
from  which  the  governing  differential  equations  rigorously  follow.   The 
equations  developed  by  the  two  methods  are  equivalent. 


The  five  diffenential  equations  from  Appendix  A  and  B  are: 
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These  five  differential  equations  and  the  other  supporting  equa- 
tions listed  in  Appendix  B  which  calculate  variables  such  as  the  mass 
transfer  flux,  humidity  ratio,  and  TEG  weight  fraction  concentration 
form  a  subroutine.   A  master  program  uses  the  forth  order  Runge-Kutta 
method  and  this  subroutine  to  produce  an  incremental  numerical  solution. 

Inputs  number  19.  The  geometry  of  the  contactor  is  specified  with 
eleven  inputs.  Initial  temperatures  take  three  inputs.  Flow  rates  and 
concentrations  each  take  two  inputs  and  a  final  input  controls  the  in- 
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crement  size  and  the  stability  of  the  solution. 

A  substantial  portion  of  the  subroutine  calculates  initial  values 
of  approximately  100  variables  including  the  fluid  properties  and  trans- 
fer parameters  discussed  in  Appendix  C. 

A  minor  portion  of  the  subrountine  is  repeatedly  called  during  the 
actual  solution.   Even  so,  approximately  three  minutes  are  required  to 
calculate  the  2500  increments  of  a  16  row  contactor  simulation.   Each  of 
the  five  state  variables  is  calculated  four  times  per  increment  while 
the  auxiliary  variables  are  calculated  only  once  per  increment. 

The  primary  outputs  of  interest  are  the  state  variable  t  and  the 
auxiliary  variables:   Y,  humidity  ratio;  P  ,  percent  TEG  weight  fraction 
concentration;  and  Z,  the  distance  in  air  flow  direction.   Cross  plots 
of  t  and  Y  produce  operating  lines  on  psychrometric  charts  while  cross 
plots  of  P  and  Z  show  the  regeneration  of  the  TEG. 


V.   RESULTS  OF  SIMULATION 

Dehumidifier  Contactor 

In  assessing  the  performance  of  the  dehumidifier  contactor,  the 
operating  line  on  a  psychrometric  chart  furnishes  the  pertinent  infor- 
mation:  the  inlet  and  exit  conditions  and  the  path  between  them.   The 
psychrometric  chart,  Figure  10a,  displays  a  parametric  study  of  dehumidi- 
fier  operating  lines  for  three  TEG  inlet  concentrations,  92%,  94%  and 
96%  and  for  three  cooling  water  temperatures,  60°,  70°  and  80°F.   All 
simulations  begin  at  92  °F  dry  bulb  and  76  °F  wet  bulb  (approximately 
50%  relative  humidity)  which  is  a  typical  outdoor  summer  temperature 
specification  for  performance  rating  of  air  conditioning  units. 

Marks  placed  on  the  operating  lines  of  Figure  10a  represent  equal 
flow  distances  of  four  rows  (6.69  inches)  except  for  the  special  mark 
which  represents  conditions  two  rows  (3.34  inches)  from  the  inlet  (top 
of  figure).   The  parameters  specified  in  the  figure  are  the  manufacturer's 
recommendations  for  a  16  row  contactor.   For  all  nine  simulations  the 
mark  indicating  8  rows  is  very  close  to  the  16  row  mark  (within  2  to  4°F 
and  3  to  4  grains) ,  and  represents  a  50%  savings  in  contactor  acquisition 
cost.  Also,  a  4  row  contactor  using  60°F  cooling  water  is  nearly  equiva- 
lent in  performance  to  a  12  row  contactor  using  70°F  cooling  water  for 
the  same  TEG  inlet  concentration.   Such  equivalent  performance  points 
can  be  the  basis  for  optimizing  cooling  water  operating  costs  and  equip- 
ment acquisition  costs  for  a  specific  application. 

The  dilution  vs.  flow  distance  of  the  TEG  solution  for  the  Simula- 
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tiotis  of  Figure  10a  is  presented  in  Figure  10b.   In  Figure  10b 
for  the  same  inlet  TEG  concentrations,  the  lowest  exit  concentrations 
correpsond  to  the  lowest  temperature  cooling  water  simulations  and  the 
drier  conditioned  air.   The  data  on  Figures  10a  and  lOh  are  from  the 
same  simulations. 

A  correlation  of  Figures  10a  and  10b  with  the  equilibrium  humectant 
values  of  Figure  4  shows  that  the  exit  conditions  of  the  16  row  contac- 
tor, operating  as  recommended  by  the  manufacturer,  are  at  equilibrium. 
For  example,  from  Figure  10b  an  exit  TEG  solution  concentration  of  90% 
corresponds  to  an  inlet  concentration  of  94%  with  60°F  cooling  water. 
On  Figure  10a  the  94%  simulation  with  60QF  cooling  water  has  exit  con- 
ditions of  60.5°F  and  0.005  H.R.  or  35  grains.   On  Figure  4  a  dry  bulb 
temperature  of  60.5°F  and  a  0.005  H.R.  corresponds  to  a  TEG  equilibrium 
humectant  value  of  90%  (the  Union  Carbide  data  base  is  used  in  the  compu- 
ter simulation),  which  is  the  exit  TEG  concentration  of  Figure  10b. 

When  psychrometric  processes  are  graphically  displayed  in  the  liter- 
ature the  path  lines  are  often  straight.   The  path  lines  for  the  liquid 
desiccant  dehumidif ication  simulation  are  not  straight  and  can  be  divided 
into  two  regions  (Refer  to  Figure  10a.)   The  first  region  is  nearly  all 
latent  cooling  (isothermal  dehumidif ication)  in  which  the  ability  of  the 
cooling  water  to  perform  sensible  cooling  (dry  bulb  temperature  change) 
is  inundated  by  the  exothermic  condensation  process.   In  the  second  reg- 
ion sensible  cooling  is  more  pronounced. 

The  first  region  physically  represents  the  flow  through  approxi- 
mately the  first  two  rows  of  a  16  row  finned  tube  contactor  or  approxi- 
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mately  3.34  inches  through  a  contactor  26.75  inches  high.   The  combin- 
ation of  a  high  rate  of  mass  transfer  at  the  inlet  (see  Figure  11)  and 
the  associated  heat  of  condensation  partially  explain  the  nearly  iso- 
thermal region.   The  exothermic  heating  due  to  mass  transfer  is  obvious 
in  Figure  12  which  shows  the  variation  of  the  liquid,  interface,  and  air 
temperatures  against  distance. 

Chemical  engineers  frequently  use  a  graph  of  the  operating  line  and 
equilibrium  curve  in  analyzing  the  absorption  of  a  substance  in  gas- 
liquid  operations.   Figure  13  shows  these  lines  for  the  present  system. 
Some  simplifications  of  the  calculation  procedure  can  be  made  if  the 
two  lines  are  straight  or  are  straight  and  parallel  [23].   However,  for 
the  LDD  system  these  lines  are  not  straight  because  of  the  temperature 
influence  on  the  equilibrium  vapor  pressure  and  hence  a  computer  simu- 
lation is  required. 

The  operating  line  of  Figure  13  represents  the  local  mole  fraction 
of  water  vapor  in  the  air  at  the  local  mole  fraction  of  water  in  the  TEG 
solution.   The  equilibrium  curve  represents  the  equilibrium  mole  fraction 
water  vapor  in  air  at  the  local  mole  fraction  water  in  the  TEG  solution. 

The  simulations  for  Figure  10a  were  iterated  to  achieve  identical 
inlet  and  exit  TEG  solution  temperatures.   This  corresponds  to  recircu- 
lating the  solution  through  the  contactor.   Thus,  the  sump  temperature 
for  each  simulation  is  slightly  less  than  the  air  outlet  temperature. 
The  greatest  change  in  liquid  temperature  is  near  the  inlet  (refer  to 
Figure  12)  so  that  the  liquid  temperature  at  the  8  row  point  of  a  16  row 
contactor  differs  little  from  the  liquid  temperature  at  the  exit.   For 
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this  reason  the  8  and  12  row  points  may  be  accepted  as  exit  conditions 
even  though  they  are  not.   Figure  14  verifies  this  point  as  the  8  row 
exit  marks  obtained  by  iterating  the  simulation  for  equal  inlet  and  exit 
liquid  temperatures  nearly  coincide  with  the  8  row  marks  of  Figure  10a, 
the  16  row  contactor. 

Figure  15  shows  the  results  of  an  air  flow  rate  which  is  75%  of  the 
air  flow  rate  of  the  Figure  14  simulations.   Slightly  drier,  1.5  grains, 
and  slightly  cooler,  1°F,  air  is  the  result  of  a  reduced  air  flow  rate 
on  the  exit  conditions. 

Specifying  a  higher  aspect  ratio  (defined  as  a  greater  face  area  to 
flow  distance  ratio  for  a  given  volume  of  contactor)  will  allow  the 
treatment  of  more  air.   Figures  10a  and  14  support  this  statement.   For 
example,  the  16  row  contactor  with  a  48"  by  127"  face  area  could  be 
arranged  as  an  8  row  contactor  with  a  96"  by  127"  face  area.   This  high- 
er aspect  ratio  contactor  will  process  twice  as  much  air  assuming  con- 
stant TEG  recirculation  rates  and  air  flow  rates  per  unit  face  area  to 
exit  conditions  within  2  to  4  °F  and  within  2  to  4  grains  of  the  16  row 
contactor  exit  conditions. 

An  indicative  result  of  an  8  row  dehumidifier  contactor  (from  Fig- 
ures 4  and  10a)  is  that  with  70°F  cooling  water  and  92%  inlet  TEG  con- 
centration the  inlet  air  at  92°F  and  111  grains  exits  at  72.5°F  and  56 
grains  or  a  dew  point  temperature  of  51°F. 

Regenerator  Contactor 

In  assessing  the  performance  of  the  regenerator  contactor  the  exit 
TEG  concentration  is  the  significant  parameter.   Figure  16a  shows  the 
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ability  of  the  regenerator  contactor  to  concentrate  90%  and  92%  TEG  so- 
lution at  temperatures  of  120°F  and  140°F.  At  constant  temperature  re- 
generation the  exit  concentration  variations  are  less  than  the  inlet 
concentration  variations.   For  example,  at  140°F  regeneration,  the  in- 
let concentration  variation  for  two  simulations  is  2%  while  the  exit 
variation  is  1/2%.   Figure  16b  shows  the  regenerator  operating  lines  on 
a  psychrometric  chart.   Both  Figures  16a  and  16b  are  for  the  same  simu- 
lations. 

The  regenerator  exit  conditions  are  nearly  at  equilibrium  as  a 
study  of  Figures  16a,  16b  and  4  will  verify.   For  example,  the  nearly 
95%  exit  concentration  of  Figure  10a  for  90%  TEG  inlet  concentration  and 
120°F  heating  water  corresponds  in  Figure  10b  to  118°F  and  0.021  H.R.  or 
147  grains.  On  Figure  4  a  dry  bulb  temperature  of  118°F  and  0.021  H.R. 
corresponds  to  an  equilibrium  humectant  value  of  about  95%  (Union  Carbide 
data  base).   The  second  equilibrium  relationship,  equation  (2.2),  for  a 
contact  temperature  of  118°F  and  95%  TEG  concentration  predicts  a  vapor 
pressure  of  24.6  mm  Hg  which  corresponds  to  0.0208  H.R. 

Y  =  0.622(— 2 )  =  0. 6220=H4)  =  °-208  <6-1> 

/OU  -  p  /Jj.4 

The  significant  result  of  the  8  row  regenerator  contactor  simulation 
is  that  TEG  concentrations  of  97%  to  97.5%  can  be  achieved  with  140°F 
heating  water,  and  95%  to  96%  TEG  concentrations  can  be  achieved  with 
120°F  heating  water. 
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VI.   CONCLUSIONS 

1.  Contactor  equipment  costs  may  be  substantially  reduced  because  half- 
sized  contactors  can  provide  95%  of  the  dehumidif ication  and  90%  of  the 
regenerated  TEG  concentration  as  compared  to  the  full  sized  contactors 
recommended  by  the  equipment  manufacturers. 

2.  Most  effective  dehumidif ication  per  unit  contactor  volume  occurs  at 
the  highest  ratio  of  face  area  to  flow  distance. 

3.  The  sensitivity  of  the  regeneration  process  to  the  regeneration 
temperature  is  minimal. 

4.  Performance  estimates  of  dehumidif ication  processes  which  continue 
to  equilibrium  can  be  prepared  from  Figure  4,  Triethylene  Glycol  Equili- 
brium Humectant  Values  on  a  Psychrometric  Chart.   First,  assume  that  the 
exit  conditions  lie  on  a  vertical  line  at  the  cooling  water  temperature. 
Then,  since  each  point  on  the  line  represents  a  humidity  ratio  and  a 
humectant  value,  and  the  sum  of  the  air  and  liquid  masses  is  constant, 
the  equilibrium  point  can  be  found. 
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APPENDIX  A 
CONTROL  VOLUME  DERIVATION  OF  RATE  RELATIONSHIPS 

The  equations  derived  in  this  appendix  are  concerned  with  the  inter- 
phase transfer  of  mass  (water  vapor)  and  energy  which  occurs  during  de- 
humidif ication  of  the  humid  air  by  the  concentrated  liquid  desiccant. 
Knowledge  of  the  equilibrium  characteristics  of  the  system  is  necessary 
for  an  understanding  of  the  mass  transfer  operation  and,  because  of 
simultaneous  heat  transfer,  enthalpy  characteristics  must  also  be  con- 
sidered. 

Operations  such  as  liquid  desiccant  dehumidif ication  quickly  become 
complex  owing  to  the  mass  transfer  in  both  the  gas  and  liquid  phases  and 
the  large  heat  effects  which  accompany  condensation  and  dilution.   Fur- 
thermore, the  need  to  remove  the  heat  of  condensation  results  in  a  non- 
adiabatic  operation. 

Figure  17  represents  the  direct  contact  of  the  gas  (dry  air  and 
water  vapor)  and  liquid  (TEG  and  water)  flowing  side  by  side  in  parallel 
flow  and  separated  by  the  gas-liquid  interface.   The  changes  in  tempera- 
ture, concentration,  etc.,  are  all  differential  over  the  section. 

Based  on  Figure  17,  the  balance  of  this  appendix  develops: 

(1)  An  overall  material  balance  for  control  volume  III; 

(2)  Energy  balances  for  control  volumes  I,  II  and  III; 

(3)  An  interface  control  volume  energy  balance. 

Next,  parameters  are  separated  into  catagories  and  finally,  a  solution 
sequence  is  presented. 
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Figure  17.   Differential  Section  of  Contactor 


The  interfacial  surface  area  in  the  section  is  dS.   For  the  finned 
tube  contactor  the  heat  transfer  area  per  unit  volume,  a^,  is  known 
from  calculations.   Because  of  the  low  surface  tension  of  the  TEG  solu- 
tion complete  wetting  of  the  heat  transfer  area  may  be  assumed  so  that 
the  mass  transfer  area  per  unit  volume,  a^,  is  equal  to  the  heat  trans- 
fer area: 

and, 

dS  =  a^Z  (A.  2) 

where  the  volume  per  unit  section  is  dZ. 

Consider  the  mass  transfer  operation  of  Figure  17  as  conducted  in 
a  steady  state  cocurrent  fashion  in  which  L  and  G  represent  the  super- 
ficial mass  flow  rate  of  the  non-diffusing  pure  TEG  and  dry  air  respec- 
tively.  Since  the  values  L  and  G  are  constant  through  the  dehumidifier, 
they  are  the  basis  for  the  material  balance. 

At  the  entrance  to  the  dehumidifier  the  gas  phase  contains  G  total 
lb/(hr)(sq  ft)  consisting  of  non-diffusing  G  lb/(hr)(sq  ft)  and  the 
diffusing  water  vapor.   The  water  vapor  concentration  is  given  by  Y  the 
mass  ratio,  or  equivalently  y  the  mass  fraction.   The  relationships 
among  these  variables  are: 

Gs  =  G(l-y)  =  YTT  (A*3) 

Y  =  zr* (A. 4) 

1  -  y 
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y  -  i  +  Y  (a. 5) 

Similarly,  the  liquid  consists  of  L  total  lb/(hr)(sq  ft)  containing 
x  mass  fraction  of  water  (or  mass  ratio  X)  and  L  lb/(hr)(sq  ft)  of  non- 
diffusing  TEG.   The  relationships  among  these  variables  are: 


Ls 

=  L< 
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-  X 
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(A. 6) 
(A. 7) 
(A. 8) 


The  quantity  of  water  vapor  in  the  gas  is  Gy  lb/(hr)(sq  ft).   In 
the  liquid  it  is  Lx  lb/(hr)(sq  ft).   Hence  the  overall  material  balance 
about  control  volume  III  of  Figure  17  for  the  water  is: 

In  =        Out 

Lx  +  Gy  =  (L+dL)(x+dx)  +  (G+dG) (y+dy)  (A. 9) 

expanding, 

Lx  +  Gy  =  Lx  +  Ldx  +  dLx  +  dLdx  +  Gy  +  Gdy  +  dGy  +  dGdy     (A. 10) 

neglecting  the  second  order  terms 

d(Lx)  =  -d(Gy)  (A. 11) 

or 

(A. 12) 
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The  mass  transfer  rate  in  lb  vapor/(hr) (sq  ft)  of  dehumidif ier  cross 
section)  is  [36] 


_  lb  water  vapor  absorbed 
(hr) (interfacial  surface) 


(A.13) 


(A. 14) 


Here  y*  is  the  equilibrium  mole  fraction  concentration  of  water  vapor 

over  the  TEG  solution  at  the  interface;  y   is  the  bulk  mole  fraction 

m 

concentration  of  water  vapor  in  the  gas  and  F  is  the  mass  transfer  co- 
efficient for  large  mass  transfer  rates  and  is  defined  by  Eq.  (C.l). 
Equation  (A. 14)  is  a  special  case  for  the  absorption  of  one  component. 

The  gas  sensible  heat  transfer  rate  across  the  interface  in  BTU/(hr) 
(sq  ft)  is  [36] 

or 

ti)aHdZ  (A. 16) 

where  C  =  heat  capacity  of  water  vapor  at  constant  pressure, 

BTU/(lb  mass) 
h  =  convection  heat  transfer  coefficient,  BTU/(hr)(°F) 

(sq  ft  interface)  and  is  defined  by  Eq.  (C.2). 
Here  h'  includes  the  effect  of  mass  transfer  on  heat  transfer  by 
accounting  for  the  sensible  heat  carried  by  the  diffusing  water  vapor. 
Radiation  has  been  neglected. 


The  sensible  heat  transfer  through  the  liquid  in  BTU/(hr)(sq  ft) 
is  [36] 

qsLaHdZ  =  VH(ti  "  CL)dZ  (A'17) 

Energy  balances  around  each  of  the  three  envelopes  of  Figure  17  can 
now  be  developed. 

Envelope  I   Energy  Balance 

In     =  Out 

GSHG  "  GS(HG  +  dHG)   +  "sgV2  +  NAaMdZ(CA(Vto)+Xo)  (A'18> 

where  the  last  term  is  the  enthalpy  of  the  transferred  vapor  relative 
to  a  reference  state  at  t  where  X  is  the  latent  heat  of  vaporization 
at  t  .   Since 


•m 


N^dZ  =  -d(Gy)  =  -d  [j^]     =  "GgdY  (A. 19) 

substitution  yields 

°  =  GSdHG  +  VH(CG-tl)dZ  "  GS(CA(Vto)+Xo)dY  (A-20) 

The  enthalpy  of  the  air-water  vapor  mixture  per  pound  of  dry  air  is 

(A. 21) 

where  C  =  heat  capacity  of  dry  air  at  constant  pressure, 
BTU/(lb)(°F) 


Define 


CS=CB+YCA  (A. 22) 

=  humid  heat  of  vapor-air  mixture,   BTU/(lb  dry  air) (°F) 


Differentiating  equation  (A. 21) 

C_dt_  +  YC.dt.  +  (C.(t  -t  )+  A  )dY 


dH^ 


B   G     AG 


Av  G  o' 


(A. 23) 


(A. 24) 


Combining  with  (A. 20)  gives  the  envelope  I  energy  balance 

-GSCSdtG  "  Vh  (W  dZ 

Envelope  II  Energy  Balance 

In  +  heat  generation  =       Out 

L  H  +Q  +NAa  dZCAT  (t   . 
S  L  XD  A  M   AL    i  o 


(A. 25) 


(A. 26) 


(A. 27) 


QC  =  ^C^L^f^2 

=  heat  removed  by  cooling  coil,  BTU/hr 
where      U  =  over-all  heat  transfer  coefficient  between  the 
liquid  and  cooling  water,  BTU/(hr)(sq  ft)(°F) 
a  =  cooling  heat  transfer  area  per  unit  volume 

sq  ft/cu  ft 
tf  =  cooling  water  temperature,  °F 

Substituting  equations  (A. 17)  and  (A. 19)  into  (A. 27)  yields 


(A. 28) 


%  -   Vil/W"  +  hL«H(ti-tL)dZ  =  LSdHL  +  QC  (A"29> 

Here 

QD  =  NaVZHD  =  -GsVY  (A-30> 

=  the  exothermic  heat  of  dilution  (mixing) ,  BTU/hr 
where  H  =  heat  of  dilution,  BTU/(lb  water  transferred) 
Substituting  equations  (A. 28)  and  (A. 30)  into  (A. 29)  yields 

"VVCai/VV^  +  hLaH(ti-tL)dZ  =  VV^VVV**        (A'31) 
The  enthalpy  of  the  TEG  -  water  solution  per  pound  of  TEG  is 

\  =  WV  +  XCAL(Vto)  (A-32) 

heat  capacity 
BTU/(lb)(°F) 

C.T  =  heat  capacity  of  liquid  water  at  constant  pressure 
BTU/(lb)(°F) 

Define 

CSL  "  CC  +  XCAL  (A-33) 

=  the  humid  heat  of  the  TEG-water  solution, 
BTU/(lb  TEG)(°F) 

Differentiating  equation  (A. 32) 

(A. 34) 


°r-  dHL  =  CSLdtL  =  CAL  (tL-Co)dX  <A'35> 

Using  equation  (A. 12) 

dHL  =  CSLdtL-^CAL  (tL-Co)dY  (A.  36) 

LS 

Combining  with  (A. 31)  gives  the  envelope  II  energy  balance 

LSCSLdtL=  (\VZ-GSCALdY)(ti-tL>-GSHDdY-DaC(tL-tf)dZ  (A-37> 

Envelope  III  Energy  Balance 

The  overall  enthalpy  balance  for  control  volume  III  is 
In  Out 

LSHL  +  °SHG  +  QD  "  LS(HL  +  dHL>  +  GS(HG  +  "V  +  QC  (A.38) 

Rearranging 


Substituting  into  equation  (A. 37) 


(A. 39) 


Uac(tL-tf)dZ  (A. 40) 

Interface  Control  Volume  Energy  Balance 
In  Out 

'iWM''  +  W2  "  VMdZ(CAL(trto))  +  «SlVZ    (A"41) 


Rearranging, 

HAaM(CA(tG-to>-CAL(ti-to)  +  V  '  WW-W'i"       (A-42> 
Solving  for  t. 


(A. A3) 


t,  -   "aWW^AlW+Wl^gV 
NAVAL  +  aH(W 

Parameter  Classification 

1.  Known  quantities  for  which  inlet  values  are  given 

2.  Constants 


U,  t 
o 

3.   Variables  with  known  functions 

V  V   V  V  V  V  y*'  ym 

By  starting  at  the  top  of  the  dehumidif ier  with  the  known  quantities 

and  solving  the  equations  at  each  incremental  step,  dZ,  the  following 
quantities  will  vary; 

V  Y,  HG>  G,  y,  Cg. 


V  V  V  '  .  v 

V  V 

and  will  be  obtained  after  each  step  and  at  the  exit  of  the  dehumidif ier . 

Solution  Sequence 

A  solution  sequence  based  on  the  differential  section  of  Figure 
17  is: 

V  =  r^-  (A. 44) 

1-y 

water  mass  ratio, 

X  =  z~-  (A. 45) 

1-x 

humid  heat, 

Cs  =  CB  +  YCA  (A.46) 

liquid  humid  heat, 

CSL  "  CC  +  XCAL  (A-47) 

gas  enthalpy, 

Hp  =  C_(t  -t  )  +  YA^  (A. 48) 

G     b   G   o       o 

liquid  enthalpy, 

"l  +  CSL(Vto)  <A-49) 

interface  temperature,  t.,  from  Eq.  (A. 43); 

equilibrium  vapor  pressure  above  liquid,  p*,  from  Eq.  (2.2); 

mole  fraction  concentration  in  the  gas  in  equilibrium  with  the  bulk 

liquid  from  Henry's  Law, 

y*  =  2-  (A. 50) 


mole  fraction  concentration  in  gas, 

ym  =  0.622  +  0.378  y  (A. 51) 

mass  transfer  flux, 

NA  "  FG  ln(fc?L)  (A"52) 

Jm 

convection  gas  heat  transfer  coefficient  corrected  for  simultaneous 

mass  transfer 

h'  _      NA°A  (A. 53) 

"  l-exp(-NACA/hG) 

This  completes  the  supporting  equations  for  the  solution  sequence.   The 
derivatives  for  the  solution  follow: 

^G  =  -Vh  (W  (A-M) 

dZ     GSCS 

dZ      b  bL 

Uac(tL-tf))  (A. 55) 


dZ    AfCAL 


1   Ua„(t-t.)  (A. 56) 


Cv  L  f 


dY    NA  (A. 57) 


(A. 58) 


dy_  =      1     dY  (A. 59) 

dZ 


dx      1       dX  (.60) 

dZ 


d(Gy) 
dZ 

GS   dZ 

d(Lx) 

-    c    ^ 

(A. 61) 


(A. 62) 


dG=  1  djGyl_   dy_  m_  (A.63) 

dZ   y  V  dZ     ^  dZ;      TM 

dL  _  1  ,d(Lx)N     dx. 

dz  "  x  (~dT-)  -  L  dZ>  =  ^^  (a-64) 


Before  simulation  of  the  contactor  with  the  solution  sequence  can 
proceed,  (1)  input  values  must  be  initialized  to  correspond  to  the  phy- 
sical system;  (2)  equations  for  fluid  properties  as  functions  of  temper- 
ature and  concentration  are  needed;  (3)  heat  and  mass  transfer  coeffici- 
ents for  the  physical  system  must  be  predicted  or  measured;  (4)  a  ra- 
tional selection  of  the  minimum  number  of  variables  necessary  to  specify 
the  state  of  the  model  must  be  made. 

Equations  (A. 54)  through  (A. 64)  above  are  not  all  independent. 
Appendix  B  describes  through  bond  graph  techniques  the  selection  of  the 
minimum  number  of  state  variables  to  represent  the  LDD  contactor.   The 
results  of  Appendix  B  indicate  that  differential  equationa  (A. 54),  (A. 55), 
(A. 56),  (A.63),  (A. 64)  above,  plus  necessary  auxiliary  equations,  will 
model  the  contactor. 


APPENDIX  B 
BOND  GRAPH  CONTACTOR  MODEL 
Introduction 

The  bond  graph  method  developed  by  Professor  H.M.  Paynter  of  MIT 
for  systematic  derivation  of  differential  equations  for  the  transport  of 
power  is  described  briefly:   Bond  graphs  are  a  powerful  analytical  tool 
as  they  give  a  physical  presentation  of  power  flows  by  maintaining  the 
topological  structure  of  the  system  while  allowing  separate  determin- 
ation of  the  sign  of  positive  power  flow  and  distinction  between  depen- 
dent and  independent  variables.   While  each  engineering  discipline  has 
its  own  calculation  methods,  bond  graphing  is  compatible  with  all  dis- 
ciplines and  is  invaluable  for  model  development  at  interfaces,  i.e. 
electrical  and  fluid.   Finally,  the  transformation  from  a  bond  graph  to 
equivalent  differential  equations  follows  an  exact  method. 

Bond  graph  theory  has  been  treated  rather  thoroughly  elsewhere  [13], 
[14],  [25],  [34],  [35].   Thus,  only  those  aspects  necessary  for  under- 
standing the  present  model  development  will  be  discussed  and  only  a 
meager  sampling  of  the  analytical  potential  of  bond  graph  theory  is  pre- 
sented. 

General  Bond  Graph  Model 

A  bond  graph  representation  of  a  control  volume  around  the  gas- 
liquid  interface  is  shown  in  Figure  18.   This  bond  graph  is  a  distributed 
resistance-capacitance  model  of  the  interface  and  can  be  replicated  n 
times  to  provide  a  model  of  any  desired  accuracy.   Dynamic  simulation  in 
time  and  space  is  possible  with  this  general  model. 
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Figure  18.   Bond  Graph  Representation  of  Coupled  Thermal  and  Mass  Transfer 
Operations  in  Liquid  Desiccant  Dehumidif ication 
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The  zeros  represent  parallel  junctions  at  which  the  sum  of  the 
heat  or  mass  flows  equal  zero.   A  distinct  temperature  or  chemical  po- 
tential is  associated  with  each  zero  junction.   Two  types  of  capaci- 
tances are  represented  on  the  graph:   thermal,  C  and  mass,  C  .   There 

t  m 

are  two  types  of  thermal  resistance:   convection  resistance  in  the 
fluids,  R  ,  and  interfacial  resistance,  R^.  .   Mass  transfer  resistance 
is  also  of  twc 

sistance,  R  . 

r 

The  dotted  line  is  a  reaction  coupling  from  the  mass  transfer  por- 
tion of  the  bond  graph  to  the  heat  transfer  portion  of  the  bond  graph 
and  represents  the  heat  of  condensation  and  dilution. 

Specific  Bond  Graph  Model 

Looking  at  a  time  invarient  solution  of  the  general  model  allows 
the  n  replications  to  be  done  with  a  dynamic  simulation  routine  if  the 
differentials  are  taken  with  respect  to  dZ.   Figure  19  shows  this  speci- 
fic simplified  model  with  the  addition  of  a  cooling  medium  with  heat 
capacity  C  and  temperature  t  . 

Bond  graphing  allows  a  determination  of  the  number  of  independent 
variables  by  counting  the  number  of  storage  elements  in  the  model.   The 
temperatures  t  ,  t   and  t   are  chosen  to  represent  the  mass  transfer 
independent  variables  because  of  their  physical  significance  and  their 
association  with  the  storage  elements.   Figure  20  is  a  word  representation 
of  Figure  19  and  the  contactor.   (Note  that  the  t's  represent  temperature 
differences. ) 
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Figure  20.   Word  Bond  Graph  of  Contactor 


Bond  Graph  Symbols 

On  Figure  19  the  bonds  with  arrows  represent  signal  flows  with  no 
feedback  along  the  same  bond.   For  example,  the  interfacial  resistance 
to  mass  transfer,  R. ,  is  dependent  on  the  interface  temperature.   Solid 
bonds  represent  the  power  flow  associated  with  heat  or  mass  transfer. 
Ones  represent  series  junctions  at  which  the  sum  of  the  efforts,  either 
temperature  or  chemical  potential,  is  zero. 

In  a  true  bond  graph  the  product  of  the  effort  variable  (above  the 
bond)  and  the  flow  variable  (below  the  bond)  is  power.   Figure  19  is  a 
pseudo  bond  graph  because  the  concentrations  are  used  as  the  mass  trans- 
fer efforts  rather  than  chemical  potential  and  because  heat  flux  is  used 
as  the  thermal  flow  rather  than  entropy.   With  a  pseudo  bond  graph  the 
reaction  couplings  must  be  scrutinized  as  standard  bond  graph  rules  do 
not  hold,  but  with  consistent  units,  pseudo  bond  graphs  follow  the  re- 
mainder of  the  rules  for  regular  bond  graphs.   The  pseudo  bond  graph  was 
selected  for  the  physical  significance  of  the  effort  and  flow  variables. [14] 


formed  before  use.   The  dotted  line  indicates  the  transformer  modulus 

is  controlled  by  t  .   Linearity  is  not  implied  by  the  bond  graph.   Two 

mass  flux  controlled  heat  sources,  S.1  and  S_,  represent  the  latent  heat 

rl      rd 

of  condensation  and  the  heat  of  dilution. 

Heat  flux  in  BTU/(hr)(sq  ft)  and  mass  flux  in  lb  water/ (hr) (sq  ft) 
are  the  flow  variables  used  in  Figure  19.   Table  1  lists  the  units  of 
the  thermal  bond  graph. 

Variable  Units 

Effort  °F 

Flow  BTU 


(hr)(sq  ft) 


Resistance  (hr) (sq  ft) (°F) 

BTU 


Capacitance  BTU 

(hr)(ft)(°F) 


Table  1.   Thermal  Bond  Graph  Units 

Dimensionally,  the  product  of  resistance  and  capacitance  is  (ft)  or 
the  same  as  the  solution  increment. 

Differential  Equation  Model 

The  differential  equations  which  result  from  a  rigorous  application 
of  bond  graph  techniques  have  been  presented  as  equations  (4.1)  through 
(4.5)  and  are  equivalent  to  equations  derived  in  Appendix  A.   The  com- 
puter simulation  uses:   (1)  differential  equations  (4.1)  through  (4.5); 
(2)  the  auxiliary  equation  for  t.,  Eq.  (4.6);  (3)  the  equilibrium  relation- 
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ship  for  vapor  pressure  over  TEG,  Eq.  (2.2);  and  (4)  the  following  se- 
quence of  auxiliary  equations: 
mass  fraction  in  liquid, 

L  "  Ls 

x=  =■  (B.l) 

L 

percent  weight  fraction  concentration  of  TEG, 

PT  =  100  (1  -  x)  (B.2) 

interface  temperature  from  Eq.  (A. 6), 

t.  =  fct<tg>VNA) 
equilibrium  vapor  pressure  over  TEG  from  Eq.  (2.2), 

p*  =  fct(t.,PT) 

Henry's  Law  for  the  mole  fraction  concentration  in  the  gas  in  equilibrium 
with  the  bulk  liquid, 

y*=^  (B.3) 

Pi- 


mass  fraction  in  gas, 

G  -  G0 


mole  fraction  concentration  in  gas, 


0.622  +  0.378y 


(B.A) 


(B.5) 


mass  transfer  flux  (F  is  derived  in  Appendix  C) , 

1  -  y* 
NA  =  FG  ln( )  (B.6) 

1  -  ym 


heat  of  dilution  taken  as  10%  of  the  heat  of  condensation  and  sensible 
heat  associated  with  the  transfered  mass, 

QD  =  NA  (100  +  (t.  -  t±))  (B.7) 

mole  fraction  concentration  in  liquid, 


xm  =  2 (B.8) 

8.334  -  7.334  x 


and  humidity  ratio, 


Y  «  j^j  (B.9) 


The  computer  model  is  now  complete. 


APPENDIX  C 
FLUID  PROPERTIES  AND  TRANSFER  PARAMETERS 

This  appendix  presents  the  essential  but  less  than  obvious  relation- 
ships.  The  development  of  an  equation  for  the  vapor  pressure  over  TEG 
from  graphical  data  is  discussed  in  Chapter  II  and  is  representative  of 
the  method  of  the  numerous  curve  fits  presented  here.   The  remaining  de- 
tails may  be  found  in  the  computer  program  listing.   Appendix  E. 

Essential  to  the  model  is  a  numerical  value  for  the  mass  transfer  co- 
efficient, F  .   This  coefficient  is  predicted  from  the  more  easily  ob- 
tained gas  convection  heat  transfer  coefficient  by  means  of  the  well 
accepted  Chilton-Colburn  analogy.   Agreement  of  the  analogy  with  newer 
data  on  mass  transfer  is  very  good  [32].   The  mass  and  gas  heat  transfer 

coefficients  are  related  by 

h       0.667 

Fp-^-S-)  (CD 

G   Cfi   Sc 

F  =  mass  transfer  coefficient,  (lb  water) /(hr) (sq  ft) 
where     C  =  heat  capacity  of  air  at  constant  pressure, 
BTU/(lb  mass) 
Pr  =  Prandtl  number 
Sc  =  Schmidt  number 

h  =  heat  transfer  coefficient  between  the  gas  and 
laminar  TEG  film,  BTU/(hr)(sq  f t) (°F) 

A  numerical  value  for  h  is  based  on  an  empirical  correlation  for  a 
representative  finned  tube  extended  surface  heat  exchanger  reported  by 
Rohsenow  [30]. 


n  to/  1,    n  r  °-681    °-333   °-2   0.113 

.     0.134  k    D  G  ,n   x      ,sN    ,s. 

hG  =  B   r  max       (Pr)     (j)    (-)  (C.2) 

r       yB 

where 

D  =  root  diameter  of  tube,  ft 
r 

k   =  thermal  conductivity  of  air,  BTU/(hr) (ft) (°F) 

G    =  mass  flow  rate  at  minimum  cross  section, 

max 

lb/(hr)(sq  ft) 
s  =  distance  between  adjacent  fins,  in 
1  =  fin  height  or  half  the  mean  fin  distance 

between  tubes,  in 
t  =  fin  thickness,  in 
y  =  viscosity  of  air  at  bulk  temperature,  lb/(hr)(ft) 

Eq.  (C.2)  shows  excellent  agreement  with  Kays  [15]  plate  finned  circu- 
lar tube  data  but  which  geometry  differs  from  the  modeled  contactor. 
A  curve  fit  developed  for  the  Prandtl  number  for  air  is 

Pr  =  0.721  exp(-196E-6  t  )  (C.3) 

To  compute  the  Schmidt  number 

Sc  =  JL_  (c.4) 

v 

values  for  the  viscosity,  density  and  diffusivity  are  needed. 

Curve  fits  led  to 

0.159 
un  =  0.0223(t  )  (C5) 


492  [18.02  y  +  28.96  (1-y)] 


(C.6) 


80 
while  the  diffusivity  is  [2] 

0.000146    (t     +  460)2'5 

Dv  = r-+Soi (c-7) 

g 

The  thickness  of  a  falling  film  of  the  water  TEG  solution  on 
the  vertical  aluminum  fin  is  [26] 

P  D  g 

where  6  =  film  thickness,  ft 

1-L  =  film  liquid  viscosity,  lb/(hr)(.ft) 
r  =  liquid  loading  per  unit  width  of  fin,  lb/(hr)(ft) 

p  =  density  of  film  liquid,  lb/(cu  ft) 

2 
g  =  acceleration  due  to  gravity,  ft/hr 

The  mass  flow  rate  per  unit  width  r  can  be  calculated  from  the  liquid 
recirculation  rate  but  values  for  y  and  p   for  TEG-water  solutions 
are  obtained  from  curve  fits 

f,      (1.487  -  0.727- In  P  ) 
yD  =  4.46  E-7(PT)lD-0:5  tl  L  (C.9) 

PD  =  54.88  +  3.73  In  PT  -  0.0262  t±  (G.10) 

The  liquid  phase  heat  transfer  coefficient  was  approximated  as 
conduction  dominated  between  the  air-liquid  interface  and  the  liquid 
mid-point  by 
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k 
\   =  -y  (Cll) 

from  the  curve  fit  for  the  TEG-water  solution  conductivity 

kD  =  0.789  -  0.14  In  PT  -  131E-6  t±  (C.12) 

The  heat  transfer  coefficient  from  the  liquid  mid-point  to  the 
aluminum  fin  was  approximated  by  the  empirical  correlation  for  a 
falling  film  [17],  [26] 

0  333 

55  j "» p°1,33c°  j  '     k4°-m  (C13) 

for 

4  r 

<   2100 

yD 

where  L  =  length  of  heat  transfer  surface,  ft 

CD  =  CSL  defined  ln  Ecl'  (A-33> 

A  curve  fit  for  the  specific  heat  of  the  TEG-water  solution  is 

C^  =  0.4977  +  0.000444  t.  +  XCAT  (C.14) 

D  1     AL 

A  cautionary  note  with  Eq.  (C.13),  [26]  suggested  a  range  of  validity 
on  L  as  0.4  to  6  ft.   Accordingly,  the  contactor  height  1.12  to  2.24 
ft  was  used  for  L.   This  gave  values  for  h  from  40  to  60  while  values 
for  h   ranged  from  150  to  200.   These  values  indicate  that  Eq.  (C.13) 
may  not  be  used  since  the  heat  transfer  coefficient  can  not  be  less 
than  that  of  pure  conduction.   Exclusive  use  of  the  conduction  coefficient 
h  is  correct.   The  best  approach  is  to  calculate  both  h.  and  h  and  use 
the  larger  for  the  entire  liquid  film.   Sole  use  of  h     produces  slightly 
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cooler  operating  lines  as  shown  by  Figure  21.   In  addition  the  typical 
overall  heat  transfer  coefficient  from  the  liquid  to  the  cooling  water 
increases  from  14  to  20  while  the  typical  overall  heat  transfer 
coefficient  from  the  air  to  the  cooling  water  increases  from  6  to  7. 
The  comparison  in  Figure  21  shows  that  the  simulation  results  are  not 
materially  affected. 

A  curve  fit  was  developed  from  tube  side  water  heat  transfer  co- 
efficient data  [17] 

D  -  2  6 
ht  =  (0.905  -  0.1863  In  (  £q7q833^) } (163  +1'731  V  Vf      (C*15) 

where 

6  =  tube  wall  thickness,  ft 
w 

Vf  =  water  velocity  through  tube,  ft/sec 

The  free  flow  area  of  Figure  8  is  slightly  reduced  by  the  TEG-water 

solution  film  thickness  so  that  the  ratio  of  free  flow  area  to  face 

area  is 

l-P(t+26)(Tsh-12Dr)  (cl6) 

Tsh 

where  P  =  number  of  fins  per  inch,  1/in 

T  .  =  horizontal  tube  spacing,  in 
sh 

The  overall  heat  transfer  coefficient  from  the  gas  to  the  cooling 

water  is  calculated  by 

1 

U  =  _i_  +  _J_  x  ^  6w  ,  aH (C.17) 

nh„   nh 
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k  =  conductivity  of  tube  wall,  BTU/ (hr) (f t) (°F) 

6  =  thickness  of  tube  wall,  ft 
w 

The  calculations  were  done  with  the  fin  efficiency  applied  to  both 
the  fin  area  and  the  tube  area. 

Applying  the  fin  efficiency  to  the  total  area  instead  of  just  the 
fin  area  accounts  for  97%  of  the  heat  transfer: 

n  '66      -  .97  (C.18) 


a    afc    .66(.93)+.07 

The  remainder  of  the  contactor  geometry  and  transfer  parameter 
calculations  is  straight  forward.   A  complete  listing  may  be  found  in 
Appendix  E. 


APPENDIX  D 
TRIETHYLENE  GLYCOL  PROPERTIES 

Triethylene  glycol  is  a  colorless,  combustible,  hygroscopic  water- 
soluble  liquid  used  as  a  chemical  intermediate,  solvent,  bactericide, 
humectant,  and  fungicide. 

Physical  Properties  [12],  [28]: 


Formula  CH  -  0  -  CH  - 
1 

CH  -   OH 

1 

CH2  -  0  -  CH2  - 

CH  -  OH 

Molecular  Weight 

150.2 

Boiling  Point  at  1  atm 

545.9  °F 

Vapor  Pressure  at  77°F 

<0.01  mm  Hg 

Specific  Gravity  at  77°F 

1.119  g/cc 

Freezing  point 

19  °F 

Pour  Point 

-73  °F 

Viscosity  at  77°F 

37.3  centipoises 

Surface  Tension  at  77 °F 

45  dynes/cm 

Specific  Heat  at  77°F 

0.53  BTU/(lb)(°F) 

Flash  Point 

320  °F 

Fire  Point 

330  °F 

Conductivity  at  77°F 

0.13  BTU/(hr)(ft)(°F) 

Desired  Desiccant  Properties   [33] 

A  liquid  desiccant  selected  for  a  dehumidif ication  system  should 
have  the  following  characteristics: 

1.  Suitable  vapor  pressure  characteristics;  should  be  capable  of  con- 
centration at  low  temperatures 

2.  Noncorrosive 

3.  Odorless 

4.  Chemical  stability;  should  not  break  down  or  vaporize  in  the  range 
of  interest 

5.  Low  viscosity  and  good  heat  transfer  characteristics;  must  not 
crystallize  or  solidify  near  the  operating  range 

6.  Nontoxic  and  nonflammable 

7.  Widely  available  and  inexpensive 

Triethylene  Glycol  Desiccant  Properties  [12],  [28] 

1.  TEG  has  a  suitable  vapor  pressure  curve  and  can  be  concentrated  at 
temperatures  as  low  as  120  °F 

2.  TEG  is  compatible  with  steel  if  slight  TEG  discoloration  due  to  iron 
contamination  is  acceptable.   Otherwise  stainless  steel,  aluminum, 
or  resin  coatings  should  be  used. 

3.  TEG  is  practically  odorless. 

4.  TEG  is  a  superior  choice  considering  both  its  hygroscopicity  (item  1) 
and  its  high  boiling  point. 

5.  TEG  has  a  low  viscosity  and  good  heat  transfer  characteristics. 

6.  TEG  is  low  in  chronic  oral  toxicity.   It  has  been  fed  in  the  drinking 


water  of  rats  over  a  two  year  period  at  a  4%  concentration  with  no 
apparent  ill  effect.   Prolonged  inhalation  of  saturated  vapors  of 
TEG  has  produced  no  ill  effects  on  humans.   Flash  and  fire  points 
determined  by  the  ASTM  Cleveland  Open  Cup  method  for  pure  TEG  are 
320  °F  and  330  °F  respectively.   TEG  presents  a  minimal  fire  hazard 
in  storage  or  handling. 
7.   TEG  is  a  byproduct  of  automotive  antifreeze.   The  cost  of  TEG  com- 
pares to  the  cost  of  antifreeze. 


APPENDIX  E 

COMPUTER  PROGRAM  LISTING 

AND 

SAMPLE  OUTPUT 
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